EPJ manuscript No. 

(will be inserted by the editor) 



Nucleon Polarizabilities: Theory 

B. Pasquini 1,a , D. Drechsel 2 • b , and M. Vanderhaeghen 2 ' 

1 Dipartimento di Fisica Nucleare e Teorica, Universita degli Studi di Pavia, 
and INFN, Sezione di Pavia, 1-27100 Pavia, Italy 

2 Institut fur Kernphysik, Johannes Gutenberg-Universitat, D-55099 Mainz, Germany 



Abstract. We review recent developments in the theoretical investiga- 
tion of the nucleon polarizabilities. We first report on the static po- 
larizabilities as measured in real Compton scattering, comparing and 
interpreting the results from various theoretical approaches. In a second 
step, we extend the discussion to the generalized polarizabilities which 
can be accessed in virtual Compton scattering, showing how the infor- 
mation encoded in these quantities can provide a spatial interpretation 
of the induced polarization densities in the nucleon. 

1 Introduction 

The polarizabilities of a composite system such as the nucleon are elementary struc- 
ture constants, just as its size and shape. They can be accessed experimentally by 
Compton scattering processes. In the case of real Compton scattering (RCS), the 
incoming real photon deforms the nucleon, and by measuring the energy and angu- 
lar distributions of the outgoing photon one can determine the induced current and 
magnetization densities. The global strength of these densities is characterized by the 
nucleon dipole and higher order (quasi-static) polarizabilities. In contrast, the virtual 
Compton scattering (VCS) process is obtained if the incident real photon is replaced 
by a virtual photon. This virtuality of the photon gives access to the generalized 
polarizabilities (GPs) and allows us to map out the spatial distribution of the polar- 
ization densities. 

Over the past years, the nucleon polarizabilities have been the subject of extensive 
experimental and theoretical studies (see, e.g., Refs. |l|2!3)4j for comprehensive re- 
views). In this work we summarize recent developments of the theoretical studies, and 
we refer to [5] for the corresponding analysis of the experimental work. In Sec. [2] we 
discuss the static polarizabilities and compare the predictions from different theoret- 
ical investigations. Several approaches to extract the polarizabilities from RCS data 
are presented in Sec. [3j The VCS process and the role of the GPs are reviewed in 
Sec. |4j and in Sec. [5] we discuss how to obtain a spatial representation of the infor- 
mation contained in the GPs. A short summary and some conclusions are given in 
Sec. El 
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2 Polarizabilities in RCS 

The physical content of the nucleon polarizabilities can be visualized best by effective 
multipole interactions for the coupling of the electric (E) and magnetic (H) fields 
of a photon with the internal structure of the nucleon. When expanding the RCS 
amplitude in the energy of the photon, the zeroth and first order terms follow from 
a low-energy theorem and can be expressed solely in terms of the charge, mass, and 
anomalous magnetic moment of the nucleon. The second order terms in the photon 
energy describe the response of the nucleon's internal structure to an electric or 
magnetic dipole field. They are given by the following effective interaction : 

tfiff = [5 <*E1 E 2 + \ PMI H 2 ] , (1) 

where the proportionality coefficients are the electric (oei) and magnetic (/?mi) scalar 
dipole polarizabilities, respectively. These global structure coefficients are propor- 
tional to the electric and magnetic dipole moments of the nucleon which are induced 
by the applied electric and magnetic fields. The polarizabilities have been measured 
extensively by use of unpolarized Compton scattering. A global fit to all modern 
low-energy proton Compton scattering data leads to the results [5] 

a p E1 = 1 2.1 ± 0.3 (stat.) =F 0.4 (syst.) ± 0.3 (mod.) , 

(3 P M1 = 1.6 ± 0.4 (stat.) ± 0.4 (syst.) ± 0.4 (mod.) . (2) 

Here and in the following the scalar polarizabilities are given in units of 10 -4 fm 3 , and 
the indicated errors denote the statistical, systematical and model-dependent errors. 
Equation [2] shows the dominance of the electric polarizability a E1 . The tiny value 
of the magnetic polarizability, fi p M1 , originates from a strong cancelation between 
the large paramagnetic contribution of the N — > A spin-flip transition and a nearly 
equally large diamagnetic contribution, mostly due to pion-loop effects. 

The internal spin structure of the nucleon appears at third order in an expansion 
of the Compton scattering amplitude. It is described by the effective interaction 



H$ = -4. 



\lE\E\ o ■ (E X E) + \~iM\M\ <T • (H X H) 
— 7M1B2 Eij (JiHj + 7B1M2 Hij (JiEj } , (3) 

which involves one derivative of the fields with regard to either time or space, E = 9<E 
and — ^(S7iEj + VjEi), respectively. The four spin (or vector) polarizabilities 
Keiei, 7M1M1) 7Afi£2, and 7eim2 describing the nucleon spin response at third order, 
can be related to a multipole expansion [7], as is reflected in the subscript notation. For 
example, 7mib2 corresponds to the excitation of the nucleon by an electric quadrupole 
(E2) field and its de-excitation by a magnetic dipole (Ml) field. Expanding the 
Compton scattering amplitude to higher orders in the energy, one obtains higher 
order polarizabilities to the respective order, e.g., the quadrupole polarizabilities at 
fourth order [715] , 

On the experimental side, much less is known about the spin polarizabilities, 
except for the forward (70) and backward (7^) spin polarizabilities of the proton, 
given by the following linear combinations of the polarizabilities in Eq. ([3]): 

70 = — 1E1E1 — TMlAfl ~ 7B1M2 — 7A/1_E2 , (4) 
7tt = — 1E1E1 + 7M1M1 — 1E1M2 + 1MXE2 ■ (5) 



The forward spin polarizability has been determined by the Gerasimov-Drell-Hearn 
sum rule experiments at MAMI and ELS A [9 10 11 . A recent analysis of these data 
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yields the value [39] 



7o = -0.90 ± 0.08 (stat.) ± 0.11 (syst.) , 



(0) 



where here and in the following all spin polarizabilities are given in units of 10~ 4 fm 4 . 
The recent experimental value for the backward spin polarizability has been obtained 
by a dispersive analysis of backward-angle Compton scattering [4] , 



This value includes both the dispersive and the large 7r°-pole contributions. In the 
analysis of Ref. [I], the latter contribution takes the value 7J ~p o1c = —46.7, which 
leads to 7 f sp = 8.0 ± 1.8. 

The effective Hamiltonians of Eqs. ([T]) and ^ describe a shift in the nucleon 
energies at second order in the electromagnetic (cm.) fields. Several lattice QCD col- 
laborations have recently implemented this fact as a tool to "measure" the nucleon 
polarizabilities by calculating the mass shifts in a constant background field and isolat- 
ing the quadratic response. In particular, it has been possible to determine the electric 
polarizability of the neutron and other neutral octet and decuplet baryons [12] as well 
as the magnetic polarizability of the proton, the neutron, and all other particles in 
the lowest baryon octet and decuplet states [T3]- In Ref. [T3], fi p M1 nas been calculated 
by use of the Wilson action in the pion mass range 0.5 < < 1 GeV and neglecting 
disconnected loop diagrams. For the smallest calculated pion mass of ~ 500 MeV, 
a value of f3 M1 — 2.36 ± 1.20 was obtained. While it is encouraging to see that the 
lattice result is in the right ball park when compared to the experimental value of 
Eq. ([2]), a more precise comparison clearly requires a dynamical fermion calculation, 
including disconnected loop diagrams and much smaller pion masses. Such small pion 
masses would then allow one to extrapolate safely to the physical pion mass within 
the framework of chiral perturbation theory. 

To determine the spin polarizabilities, Eq. (J3J can likewise be used to calculate 
energy shifts of a polarized nucleon in an external field. As an example, consider a 
nucleon polarized along the z-axis and apply a magnetic field rotating with angular 
frequency u in the xy plane, 



where e, stands for the unit vector in the direction i — x,y and B is the magnitude 
of the field. Such a field leads to an energy shift AE = ^2ttjmimi^Bq if the nucleon 
spin is oriented along the positive (— ) or negative (+) z-axis. The split between 
the two levels is then directly proportional to the magnetic dipole spin polarizability 
7jwimi ■ R has been shown in Ref. [O] that allowing for background fields with suitable 
variations in space and time, lattice QCD should be able to calculate all six dipole 
polarizabilities. In particular, calculations are in progress to determine the electric 
polarizability of a charged particle such as the proton as well as the four proton spin 
polarizabilities of Eq. ([3| TS] . 

A microscopic understanding of the nucleon's polarizabilities requires to quantify 
the interplay between resonance contributions, e.g., the N — > A transition, and long 
range pion-cloud effects. Systematic studies of such pion-cloud effects became possible 
with the development of chiral perturbation theory (ChPT), by an expansion of the 
lagrangian in the external momenta and the pion mass ( "p-expansion" ) . The first 
such calculation at leading order, 0(p 3 ), yielded the following values for the proton 
scalar polarizabilities [IB] : 



7tt = 



38.7 ± 1.8 (stat. + syst.) . 



(7) 



H = Bq [cos(wi) e x + sin(cji) e y ] , 



(8) 



a 



P = 

El 



mi 




= 12.2 



(9) 



Ml — 



96irf%m. 
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Table 1. Theoretical predictions for the dispersive contribution to spin polarizabilities of 
the proton: to 0(p 3 ) in HBChPT 17:, to 0(p 4 ) in HBChPT [18], to C(e 3 ) in the small scale 
expansion [19], in the Lorentz covariant (LC) ChPT to 0(p 4 ) [21 j . in the chiral approach with 
unitarity and causality (GLP) of Ref. [20] , in fixed-t dispersion relations analysis (HDPV) 
of Ref. [E] and Ref. [7] (BGLMN), and in hyperbolic DRs of Ref. [I] (HYP. DRs). 





o( P 3 ) 




<D(e 3 ) 


LC 


GLP 


HDPV 


BGLMN 


HYP. DRs 




-5.7 


-1.4 


-5.4 


-2.8 


-3.7 


-4.3 


-3.4 


-3.8 


7A/1M1 


-1.1 


3.3 


1.4 


-3.1 


2.5 


2.9 


2.7 


2.9 


7E1M2 


1.1 


0.2 


1.0 


0.8 


1.2 


-0.01 


0.3 


0.5 


7M1E2 


1.1 


1.8 


1.0 


0.3 


1.2 


2.1 


1.9 


1.6 


70 


4.6 


-3.9 


2.0 


4.8 


-1.2 


-0.7 


-1.5 


-1.1 


7^ 


4.6 


6.3 


6.8 


-0.8 


6.1 


9.3 


7.8 


7.8 



where a em = 1/137, qa — 1-27, and = 92.4 MeV. This result is in remarkable 
agreement with the experimental result of Eq. |2| . It also illustrates that these quan- 
tities diverge in the chiral limit, which is a challenge for the lattice QCD calculations. 
Conversely, the chiral expansion converges well in the "small m^" regime, and there- 
fore ChPT can complement the lattice calculations by extrapolation to the physical 
pion mass. The ChPT work was extended to 0(p 4 ) within the heavy-baryon expan- 
sion [22], yielding a E1 = 10.5 ± 2.0 and fi v M1 = 3.5 ± 3.6. The error bars for these 
values indicate that several low-energy constants had to be determined by resonance 
saturation, e.g., by putting in phenomenological information about resonance and 
vector mesons. Since the Z\(1232) lies close in energy, it has been proposed to include 
the resonance dynamically. This leads to an additional expansion parameter, the N A 
mass splitting (e expansion). Unfortunately, at C(e 3 ) the "dynamical" A increases the 
polarizabilities to values far above the data, a E1 = 16.4 and ff MX = 9.1. This can be 
changed by introducing large low-energy constants within a higher-order calculation, 
however, at expense of losing the predictive power |23j . More recently Pascalutsa and 
Phillips [24] proposed an alternative treatment in which the power counting changes 
between the threshold and the resonance regions, depending on which scales are en- 
hanced. At low energies they introduce the ratio S, with S — (Ma ~ M^)/Mn as 
new expansion parameter. In this scheme, the values for the scalar polarizabilities 
at N 3 LO are a E1 = 10.8 ± 0.7 and fi v M1 = 4.0 ± 0.7, where the error bar is an esti- 
mate of higher-order contributions [2 5) , and show a much improved agreement with 
experiment. 

Results within different theoretical approaches for the dispersive contribution to 
the spin polarizabilities are collected in Table [I] The predictions from unsubtracted 
fixed-i dispersion relations (DRs) are based on different inputs for the pion pho- 
toproduction amplitudes, while the results from hyperbolic DRs are obtained from 
dispersion integrals which run along a path at fixed angle. The agreement between the 
different DR results in Table [T] is quite satisfactory in all cases, and the spread among 
the predicted values can be seen as the best possible error estimate of such calcula- 
tions to date. The ChPT predictions in the heavy baryon, small-scale expansion and 
Lorentz-covariant approach disagree in some cases, both among each other and with 
the DR results. In particular, issues about the convergence of the chiral expansion for 
the spin polarizabilities deserve further studies. In Table [T] we also quote the results 
from a recent approach based on the chiral Lagrangian at order 0(p 3 ) involving pion, 
nucleon and photon fields, and without isobar fields as explicit degrees of freedom 
(see column labeled "GLP") [2D]. The physics of the isobar resonances is included 
by infinite summation of higher order counter terms in the chiral Lagrangian, taking 
into account the constraints from unitarity and causality. In this framework, the au- 
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thors of Ref. |26| were able to give an unified description of photon- and pion-nucleon 

scattering data up to and beyond the isobar region, with the same set of parameters 
adjusted to pion-nucleon scattering data. For jeim2 and 7mie2j the ChPT results 
at 0(p 3 ) are recovered within small numerical deviations (compare columns 0(p 3 ) 
and GLP in Table [T]), which are due to different values used for the parameters used 
in the calculations. In contrast, the ChPT and GLP results for jeiei and Ymimi 
are at variance because of large rescattering effects, and the GLP values are in much 
better agreement with the predictions from DRs. The GLP scheme shows distinct 
indications of an improved convergence pattern as compared to ChPT, as was also 
observed for neutral-pion photoproduction |26j . However, these expectations for the 
spin polarizabilities need to be confirmed by explicit calculations. 



3 Extraction of RCS polarizabilities 

The extraction of polarizabilities from RCS data has been performed by three tech- 
niques. The first one is a low-energy expansion (LEX) of the Compton cross sec- 
tions. Unfortunately this procedure is only applicable at photon energies well below 
100 MeV, which makes a precise extraction a rather challenging task. The sensitiv- 
ity to the polarizabilities is increased by measuring Compton scattering observables 
around pion threshold and into the Z\(1232) resonance region. A second formalism 
which has been successfully applied to Compton data is based on DRs, in particular 
unsubtracted [27] and subtracted [55] fixed-f DRs were used to extract the scalar po- 
larizabilities of Eq. ( 2]) . A third approach has been developed within the framework 
of chiral effective field theories [23] 24 25 29 30]. For energies below and around pion 
threshold the full Compton scattering process has been calculated to 0(p 4 ), allowing 
for an independent extraction of the polarizabilities from Compton scattering data. 
The thus obtained values for and fi v M1 in the work of |23|29j are nicely compatible 
with the results given by Eq. pj). Recently, the "<5-counting" scheme has been applied 
to extract the scalar polarizabilities using a subset of low-energy data up to E 1 = 149 
MeV |25) , and a preliminary study in a larger energy range has been presented in [31] . 
Both these analyses tend to favor a larger value for fi v M1 than given by Eq. ([2| . These 
findings call for new dedicated measurements below pion threshold. In particular, it 
has been suggested that the sensitivity to the scalar polarizabilities can be enhanced 
by combination of unpolarized cross sections and data obtained with linearly polar- 
ized photons in particular kinematics. Proposals for such measurements below pion 
threshold have been presented at the HI7S facility at Duke 32J. 

Because a reliable data analysis is based on DRs, we recall some pertinent features 
of this technique in the following. The T matrix of RCS can be expressed by 6 inde- 
pendent structure functions A;(f, t) 27J. These functions depend on the variables v 
and t, which are related to the initial (E 7 ) and final (E' ) photon laboratory energies, 
and to the scattering angle fliab by t — —AE^E' sin 2 (#i a b/2) and v — E 1 + i/(4Mjv). 
The invariant amplitudes Ai are free of kinematical singularities and constraints, they 
also obey crossing symmetry and gauge invariance. Assuming further analyticity and 
an appropriate high-energy behavior, these amplitudes fulfill unsubtracted DRs at 
fixed i, 



where A?° c is the nucleon pole term, vq is the pion production threshold and V 
denotes the principal value integral. The latter can be calculated if the absorptive 
part of the amplitude, Im Ai, is known to a sufficient accuracy. Because of the en- 
ergy weighting, the pion production near threshold and the mesonic decay of the 




(10) 



6 



Will be inserted by the editor 



low-lying resonances yield the biggest contributions to the integral. With the existing 
information on these processes and some reasonable assumptions on the lesser known 
higher part of the spectrum, the integrand can be constructed up to cm energies 
W « 2 GeV. However, Regge theory predicts that the amplitudes A\ and A 2 do not 
drop sufficiently fast to warrant a convergence of the integral. This behavior is mainly 
due to fixed poles in the t channel. In particular the i-channel exchange of pions and 
a mesons leads to the bad convergence for A 2 and A\ , respectively. The latter meson 
has a mass of about 600 MeV and a very large width, it models correlations in the 
two-pion channel with spin and isospin zero and positive parity. In order to obtain 
useful results for these two amplitudes, [27J proposed to close the contour integral in 
the complex plane by a semi-circle of finite radius v max , and to replace the contribu- 
tion from the semi-circle by a number of energy independent poles in the t channel. 
This procedure is relatively safe for A 2 because the ir° pole or triangle anomaly is 
well established by both experiment and theory. However, it introduces a considerable 
model-dependence for A\. 

In order to avoid the convergence problem and the phenomenology necessary to 
determine the asymptotic contributions, it was suggested to subtract the DRs at 
v = [28j . This subtraction improves the convergen ce b y two additional powers of 
v' in the denominator of the dispersion integrals, Eq. (10). The subtraction functions 
Ai(v = 0,i) can be obtained from subtracted DRs in t with the imaginary part of 
the amplitude 77 — > tttt — > NN as input. In a first step, a unitarized amplitude for 
the 77 — > tttt subprocess is constructed from the available experimental data. This 
information is then combined with the tttt — ¥ NN amplitudes determined by analytic 
continuation of ttN scattering amplitudes |33j . Once the t dependence of the subtrac- 
tion functions Ai(0,t) is known, the subtraction constants etj = A;(0, 0) have to be 
fixed. The 6 subtraction constants a\ to a$ are given by linear combinations of the 
scalar and spin polarizabilities and can be used to fit the RCS data. In the analysis of 
unpolarized cross sections it is sufficient to fit a\ and a 2 , or equivalently (olei ~ Pmi) 
and 7^ to the data. The remaining 4 subtraction constants can be calculated through 
an unsubtracted dispersion integral. 

In Fig. [l] we show the differential cross sections measured at several laborato- 
ries, as a function of the lab photon energy and at different scattering angles. The 
solid and dashed lines correspond to the results from unsubtracted and subtracted 
fixed-i DRs, respectively, as obtained for fixed values of OeI) Pmi and j„. Except 
for extreme backward scattering, both approaches lead to nearly identical results. A 
fit to the data within the unsubtracted DR formalism yields the values in Eq. ^ 
for the scalar polarizabilities, and compatible results, within error bars, are also ob- 
tained from subtracted DRs [T]. Because both formalisms provide an independent 
cross-check for the extraction of the polarizabilities, we can conclude that the present 
analysis of the low-energy data is well under control. We also note that very similar 
results for the cross sections below pion threshold were recently obtained within the 
unitary approach of Ref. [2U] . 

Because the unpolarized cross section is mainly sensitive to the scalar polarizabil- 
ities and to the backward spin polarizability, double polarization experiments will 
be required to get information on the individual spin polarizabilities. First measure- 
ments of these observables have recently started at MAMI The proposal is to 
measure three different observables, i) the photon-beam asymmetry S3 with linearly 
polarized photons, ii) the beam-target asymmetry E 2z with circularly polarized pho- 
tons and a longitudinally polarized target, and Hi) the beam-target asymmetry E 2x 
with circularly polarized photons and a transversely polarized target. The simulated 
data points with the expected error bars for 100 hours running time are shown in 
Fig [2] in comparison with the predictions from subtracted DRs obtained for different 
values of jeiei and Jmimi [SS]- From a preliminary analysis in Ref. [33], the pro- 
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Fig. 1. Differential cross section for Compton scattering off the proton as a function of 
the lab photon energy _B 7 and at different scattering angles. Solid lines: results from fixed- 
t subtracted DRs, dashed lines: fixed-f unsubtracted DRs. All results are shown for fixed 
values of agi + Pmi = f3.8, asi — Pmi = 10, and j v — —37. The data are from Ref. [6] 
(full circles), Ref. (diamonds), Ref. [3Z] (triangles), and Ref. fHH] (open circles). 



jected errors of the extracted spin polarizabilities are expected to be smaller than the 
spread of values presented in Table [T] Therefore, these experiments hold the promise 
to discriminate between the different theoretical predictions. 

However, all the mentioned approaches for the extraction of the polarizabilities will 
necessarily contain some model dependence through parameters describing the high- 
energy regime. Such parameters are the low-energy constants appearing in ChPT, and 
in DRs the extrapolation of the photoproduction data to regions not covered by the 
experiment. Any additional confirmation of the theoretical framework by new alter- 
native theoretical developments or independent experimental information is therefore 
most welcome. 

To this aim, the DR calculation for the double polarization observables to be mea- 
sured at MAMI, was recently compared with the unitary and causal approach based 
on the chiral Lagrangian of Refs. [20 26 . Although the two approaches are quite dif- 
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Fig. 2. The beam-target asymmetries Ei x (left), ^2z (center), and the beam asymmetry S3 
(right), as function of the cm scattering angle at fixed photon energy = 240 MeV. The 
different lines show the predictions of subtracted DRs using the experimental values for asi, 
Pmi, 70, and 7^ as given by Eqs. pj, |6| and Q, while the remaining polarizabilities are 
free parameters. Left panel: results for fixed tmimi = 2.9 and values of ■Jeiei as indicated. 
Central and right panels: results for fixed 7bibi = —4.3 and values of 7mimi as indicated. 
The plots are taken from the MAMI proposal [34], with the points representing simulated 
data with the expected error bars for 100 hours running time. 



ferent in construction, they agree well below the pion-production threshold and start 
deviating slightly when the energy increases. The small difference between the two 
calculations at higher energies can be used to estimate the theoretical uncertainty in 
calculations of Compton scattering observables. 

A second cross check was performed in Ref. j39. , by comparing the "experimental" 
spin-dependent forward Compton scattering amplitude (as constructed from the avail- 
able data on the helicity-dependent absorption cross sections ) with the DR calcu- 
lation obtained from different phenomenological inputs for the pion-photoproduction 
multipoles. The spin-dependent forward Compton scattering amplitude can be para- 
metrized as 

where the leading term is due to intermediate nucleon states (Born terms) and the 
higher-order contribution is described by the energy-dependent (dynamic) forward 
spin polarizability (FSP) 7g yn (f). This polarizability is a complex quantity which 
obeys the following DRs: 

Re[7 ° M] = v L ,v 2 - ^ ' ( } 

Im[7o dy >)] - . (13) 

For v < vq the imaginary part vanishes, and the dynamical FSP has a LEX with the 
leading-order contribution given by the forward polarizability 70. The "experimen- 
tal" dynamical FSP is compared in Fig. [3] with the predictions based on different phe- 
nomenological parametrizations of the pion photoproduction multipoles [40141142143] . 
The onset of S'-wave pion production at _E 7 = v Q leads to a strong cusp effect in the 
real part. The rapid increase of the dynamic FSP from its static value —0.90 • 10 -4 
fm 4 at E 1 = to about 6 • 10~ 4 fm 4 at pion threshold clearly shows the necessity to 
analyze Compton scattering within the framework of dispersion analysis. In particu- 
lar, such an approach is prerequisite in order to determine the 4 spin polarizabilities, 
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Fig. 3. The real (upper panel) and imaginary (lower panel) parts of the dynamic forward 
spin polarizability 7o yn as function of the photon lab energy E~ t . The light grey (yellow) and 
dark grey (green) bands show the statistical and systematic uncertainties. The predictions are 
based on the pion photoproduction multipoles of Hanstein el al. [40] (solid line), SAID [41] 
(dotted line), MAID07 g2] (dashed line), and DMT g3] (dash-dotted line). 

which yield significant contributions to the cross section only for photon energies 
above 80 MeV where the LEX does not apply. Except for the minimum of Re[7Q Yn ] 
near E y — 0.23 GeV, the experimental analysis is in very good agreement with the 
predictions of the presented models. The deviation in the minimum is not too surpris- 
ing, because this comes about by a delicate balance between the positive contribution 
from the S-wave background and the negative contribution from the low-energy tail 
of the Z\(1232) resonance. 

4 VCS and Generalized Polarizabilities 

The VCS process on the proton is accessed through the ep — > epj reaction. In this 
process, the final photon can be emitted either by the proton, which is referred to 
as the fully virtual Compton scattering (FVCS) process, or by the lepton, which is 
referred to as the Bethe-Heitler (BH) process. The amplitude T ee 7 of the ep — > ep"f 
reaction is given by the coherent sum of the BH and the FVCS process : T ee 7 ~ 
T BH +T FVCS . The BH amplitude T BH is exactly calculable from QED if one knows the 
nucleon e.m. form factors. The FVCS amplitude T FVCS contains, in the one-photon 
exchange approximation, the VCS subprocess 7*p —¥ jp. The VCS amplitude can be 
further decomposed into a Born and a non-Born contribution. In the Born process, 
the virtual photon is absorbed on a nucleon and the intermediate state remains a 
nucleon, whereas the non-Born process contains all nucleon excitations and meson- 
loop contributions. The separation between Born and non-Born parts is performed 
as described in Ref. [44145] . to which we refer for the details. The behavior of the 
non-Born VCS contribution at low outgoing-photon energy E' but at arbitrary three- 
momentum of the virtual photon, can be parametrized by 6 generalized polarizabilities 
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(GPs), which depend on the virtuality Q 2 transferred by the virtual photon. The 
GPs are denoted by P<- M ' L '< M L ~> S (Q 2 ) [44146147] . In this notation, M (M') refers 
to the electric (E), magnetic (M) or longitudinal (L) nature of the initial (final) 
photon, L (Z/) represents the angular momentum of the initial (final) photon, and 
S differentiates between the spin-flip (S — 1) and non spin-flip (S = 0) character 
of the transition at the nucleon side. Assuming that the emitted real photons have 
low energies, we may use the dipole approximation (JJ = 1). For a dipole transition 
in the final state, angular momentum and parity conservation lead to 10 GPs [H]. 
Furthermore, it has been shown [46] that nucleon crossing combined with charge 
conjugation symmetry of the VCS amplitudes provides 4 additional constraints among 
the 10 GPS, which leaves only 6 independent GPs. 

In the limit of E' —> 0, the VCS experiments can be analyzed in terms of LEXs, 
as proposed in Ref. |44j . These expansions are based on a low-energy theorem stating 
that the radiative amplitude for point-like particles diverges like 1/E' for EL — > 0, 
whereas the dispersive amplitude vanishes like E' in that limit. The interference 
between the 1 / E^ contribution of BH plus Born amplitudes and the leading term of 
the non-Born amplitude, can be expressed by 3 structure functions depending on the 
GPs, 

P LL = -2V6M N G E P (LhL1)0 , 

P TT = -3GmC(p (M1 ' M1)1 - V2q p( L1 ' M2 ^) , (14) 

V I H 2. qo 

with Ge and Gm the electric and magnetic nucleon form factors. 
However, since the sensitivity of the VCS cross sections to the GPs grows with the 
photon energy, it is advantageous to go to higher photon energies, provided one can 
keep the theoretical uncertainties under control when approaching and crossing the 
pion threshold. The situation can be compared to RCS as described in Sec. |2j where 
it was shown that a DR formalism is prerequisite to extract the polarizabilities in the 
energy region above pion threshold where the observables are generally more sensitive 
to the GPs. In order to set up DRs for VCS, we describe the VCS tensor in terms of 12 
independent amplitudes (i = 1, . . . , 12) which are free of kinematical singularities 
and constraints and even under crossing. These amplitudes depend on the 3 variables 

v, t and Q 2 , with u = E^ + (t - Q 2 )/(4M N ) , and t = 2E\ (cos # lab ^E 2 + Q 2 - 

.Ely) — Q 2 . Assuming further an appropriate analytic and high-energy behavior, these 
amplitudes fulfill unsubtracted DRs in the variable v and at fixed t and Q 2 , 

ReF," B (0V,t) = fT'°(<3 2 ,M) - f?«V,t) 

+ i v r d ^j^4A, (15) 

where the pole amplitudes F? ole are obtained from the Born amplitudes at the pole 
position, that is, with all numerators evaluated at the pole. Furthermore, Im Fi are 
the discontinuities across the s-channel cuts, starting at the pion-production thresh- 
old, which can be calculated from the absorption cross sections due to 7riV, wwN, 
and heavier hadronic states. As long as we are interested in the energy region up to 
the Z\(1232), we may saturate the s-channel dispersion integral by the ttN contri- 
bution, choosing f max w 1.5 GeV as upper limit of integration and using empirical 
information from pion photo- and electroproduction as encoded in the MAID2007 
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Fig. 4. The structure functions describing unpolarized VCS on a proton compared with 
the data from MAMI (circles |49| . squares |50|). MIT-Bates (up triangles |51|). and JLab 
(stars [52]). The RCS data [6] are shown by the (black) down triangles (slightly displace d in 
Q 2 ). All the lines are shown for e = 0.645 and based on the parameterizations of Eqs. (161 



and (17 1 for the proton GPs. The lines in the upper panel correspond to different sets of 
parameters, dotted (black): A a = 0.7 GeV, C a = (GP I), solid (red): A a = 0.7 GeV, C a = 
-150 GeV~ 7 (GP II). The solid (red) line in the lower panel is obtained for Ag = 0.5 GeV. 
The dashed (blue) lines in both panels show the contributions of the spin GPs. 



parametrization [32]. However, Eq. (15) is only valid if the amplitudes drop fast 
enough such that the integrals converge. The high-energy behavior of the amplitudes 
Fi was investigated by [15] in the Regge limit [y — > oo, t and Q 2 fixed), where it 
was found that the dispersion integrals diverge for two amplitudes, F\ and F$. The 
asymptotic contributions to these amplitudes are described by the i-channel exchange 
of a and ir° mesons, respectively, and other effects such as many-pion and heavier 
intermediate states in the s-channel integral. All these effects beyond the dispersive 
irN contribution are parametrized in terms of asymptotic contributions. In particular, 
the asymptotic part of the magnetic GP is described by a dipole function 

P^ M1)0 {Q 2 ) = ^ M/1)0 (o)/(i + g 2 /4) 2 - (is) 

To describe the electric GP, we allow for an asymptotic part consisting of a sum of a 
dipole and a Gaussian, in the same vein as the parametrization proposed in [53] for 
the nucleon form factors, 

Pisy hL1) °(Q 2 ) = P£y' L1) °(0)/(l + Q 2 /A 2 a ) 2 + C a Q 4 e-W 2 - - 15 )/ - 15 . (17) 

The values at the real-photon point are fixed by the difference between the dispersive 
7rAf contribution and the empirical information for the proton electric and magnetic 
polarizabilities as obtained from RCS experiments [B], which leads to Pasy' L1 ^° (0) = 
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Fig. 5. The spin-flip GPs without the 7r°-pole contribution. The lines shown correspond 
to the predictions from several calculations. Solid: dispersive ttN contribution [3Sj, short- 
dashed: C(p 3 ) HBChPT [541551 . long-dashed: 0(p A ) HBChPT [561571 . dashed-dotted: linear 
a model [55], dotted: non-relativistic CQM |59) . For visibility, the tiny CQM results for 

p(Ll,Ll)l and p(Ml,L2)l are multiplied by a f actor 10 Q. 

-14.37 GeV -3 , and Pis"' M1) °(0) = 21.82 GeV -3 . The remaining three parameters 
Ap,A a , and C a describe the Q 2 dependence of the asymptotic parts of the spin- 
independent GPs and can be determined by a fit to the available VCS data. In the 
upper panel of Fig. |4j we show the comparison with the experimentally measured un- 
polarized structure functions Pll~ Ptt / £ (Plt), which is proportional to the electric 
(magnetic) GPs respectively, up to a small spin-GP contribution (dashed lines). The 
magnetic GP is essentially given by Plt shown in the lower panel of Fig. [4] The figure 
shows that all the data can be well described by Ap = 0.5 GeV. For the electric GP, 
a fit to the MIT-Bates and JLab data is obtained for C a — 0, and A a — 0.7 GeV 
(denoted by parameterization GP I). However, this does not describe the MAMI data 
at intermediate Q 2 , which require an additional structure, parameterized through the 



Gaussian term in Eq. (17 1. A good description of all available data is obtained for 
A a = 0.7 GeV, and C a — —150 GeV~ 7 (denoted by parameterization GP II). 
The GPs were also calculated in HBChPT at C(p 3 ) p|55] and at 0(p 4 ) p|57] . Al- 
though the HBChPT and DR results agree qualitatively in the description of the unpo- 
larized VCS structure functions, there are large differences in the spin sector. Figure[5] 
shows the predictions for the spin GPs for several different approaches. The compar- 
ison demonstrates that a satisfying theoretical description of the spin-flip GPs over a 
large range in Q 2 is still a challenging task. This calls for VCS experiments which are 
sensitive to the spin GPs. Two types of such experiments can be envisaged: (i) unpo- 
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larized VCS with variation of the transverse photon polarization e in order to separate 
the response functions Pll an d Ptt and (it) double-polarization experiments to ac- 
cess new structure functions directly related to the spin-flip GPs [HD]. A first test 
experiment for double-polarization observables was performed at MAMI [61J . Fur- 
thermore, these structure coefficients are also prominent input in calculating hadronic 
corrections to the Lamb shift and the hyperfine splitting in hydrogen [55] . 



5 Spatial representation of GPs 

In the last years, generalized parton distributions 64 65 and form factors [66167] have 
been discussed as a tool to access the distributions of partons in the transverse plane, 
and first calculations of these spatial distributions have been performed in lattice 
QCD [5H] and hadronic models [55] . Most recently, it has also been shown how the 
concept of GPs can be used to describe the spatial deformation of the charge and 
magnetization densities if the nucleon is exposed to an external quasi-static electric 
field [70171] . 

In order to arrive at a spatial representation of the information contained in the GPs, 
we consider the VCS process in a light-front frame [7T] . The VCS tensor is defined as 

H^ = -i J d 4 xe- i "- x {p , ,X , N \T[J^(x),J v (0)]\p,X N ), (18) 

with \n and p (X' N and p') the helicity and four-momentum of the initial (final) 
nucleon. We indicate the (large) light- front + component by P + = (p + +p' + )/2 
(defining a = a ±a 3 ), and choose the symmetric frame by requiring that A = p' —p 
is purely transversal, i.e. A + = 0. In the light-front frame, the + component of 



the current in ( 18 ) is a positive-definite operator for each quark flavor, allowing 
for a light-front charge density interpretation. In the low-energy limit v — > 0, the 
VCS light-front helicity amplitude H(X' J , X' N , X N ) = e'* (A^) H +u can be described 
in terms of GPs. In particular, we consider the transverse polarization component of 
the outgoing photon corresponding to an electric field E ~ iE' e' x . Any system of 
charges will respond to such an applied electric field through an induced polarization 
Po, which will be forced to align with the external field such as to minimize the energy 
— E • Pq. The linear response in the outgoing-photon energy of the helicity averaged 
VCS amplitude therefore allows one to define an induced polarization Po as 

ei*(A^) -P = ( 2 p+^ q^ h {^ ( ,X N ,X N ) \ v=Q , (19) 

where r = Q 2 /(4M^). By Fourier transforming Po, we obtain the spatial distortion 
of the charge density in an unpolarized nucleon as induced by the external electric 
field, 

Po (b) = b J°° ^- QJ X (b Q)A(Q 2 ), (20) 
where b is the transverse position, b = |b|, and b = h/b and A(Q 2 ) can be expressed 



in terms of the scalar and spin GPs. The dipole pattern described by Eq. (20) is 
shown in Fig.|6]for the parameterizations GP I and CP II in the left and right panels, 
respectively. The comparison shows that the enhancement at intermediate Q 2 in GP 



II (second term in Eq. (17)) gives rise to a much larger transverse extension of the 
induced polarization cloucL Forthcoming VCS experiments planned at MAMI [75] are 
conceived to study the Q 2 dependence of the GPs in more detail. In this way, the 
experiment is expected to verify or disprove the large distance structure predicted by 
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Fig. 6. The induced polarization Pq in a proton of definite light-cone helicity if exposed to 
an e.m. field with photon polarization along the x-axis, as indicated. The left (right) panel 
is obtained for GP I (GP II), see caption of Fig. [4] The light (dark) regions correspond to 
the largest (smallest) values. 



model GP II. 

Analogously, we can define the linear response to an external quasi-static e.m. 
dipole field if the nucleon is in an eigenstate of transverse spin, Sj_ = cos <ps e x + 
sin</>se y , with (ps the angle indicating the spin vector direction. The dependence of 
the induced polarization on the transverse position is given by 



P T (b)=Po(b)-6(S ± xe,)-6 



Jo ( 27r ) 



dQ 

(27T) 



QJ 2 {bQ)B 



MbQ)C + J -^B 
bQ 



(21) 



displaying dipole, quadrupole and monopole patterns with weights B, and C given 
by combinations of scalar and spin GPs. 

In Fig.[7]we show the spatial distributions in the induced polarization for a proton 
of transverse spin (chosen along the x-axis) for parameterization GP II. The compo- 
Pq displays a quadrupole pattern with pronounced strength around 0.5 fm 



ncnt Pf - 



due to the electric GP, whereas the component — Pq shows an additional monopole 
pattern, which is dominated by the 7r°-pole contribution. 



6 Conclusions 

In this contribution we have discussed real and virtual Compton scattering off the 
proton as powerful tools to extract information on the proton polarizabilities. In the 
case of RCS, the static scalar polarizabilities are now known with relatively small error 
bars, whereas our knowledge on the spin-dependent sector is as yet incomplete. Only 
two combinations of the spin polarizabilities have been determined as yet. However, 
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Fig. 7. The nduced polarization density in a proton, with spin S oriented along the x-axis, 
when submitted to an e.m. field. The left (right) panels is for P^ — Pq (P^ — Pq), and 
correspond with photon polarization along the x-axis (j/-axis) as indicated. The light (dark) 
regions correspond to the largest (smallest) values using parameterization GP II, see caption 
of Fig. [4] 



the forthcoming double-polarization experiments at MAMI will allow us to disentangle 
all 6 spin-polarizabilities individually. The spin polarizabilities of the nucleon are 
fundamental structure constants of the nucleon, just as shape and size of this strongly 
interacting many-body system. Therefore, the ongoing experimental programs are 
both timely and important for our understanding of the nucleon structure. These 
activities will provide stringent precision tests for both the existing predictions of 
effective field theories and future results expected from the lattice-gauge community 
for the polarizability of the nucleon. 

Moreover, the generalized polarizabilities from VCS processes allow us to map 
out the spatial dependence of the induced polarizations in an external e.m. field. The 
existing field theoretical formalism to extract light-front densities from nucleon form 
factor data has been recently extended to the deformations of quark charge densities 
under the influence of an applied e.m. field. It has been shown that the available 
proton electric GP data yield a pronounced structure in the induced polarization 
density at large transverse distances of 0.5 — 1 fm. At Q 2 values smaller than 0.1 
GeV 2 , chiral effective field theory was found to describe the VCS data quite well, 
which highlights the role of pions in the nucleon structure. Such description can 
however not be applied at intermediate and large Q 2 values. This transition region 
is dominated by nucleon resonance structure, which can be described by dispersion 
relations. Forthcoming VCS precision experiments at MAMI in this intermediate Q 2 
region will be capable to determine the spin polarizabilities in more detail and, in this 
way, help to unravel the the distribution of quark charge and magnetization currents 
in the nucleon. 
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